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Abstract: This is a short review of the physical ideas, algorithm for cal- 
culations of the phase transition pressure and some results of a semi-classical 
theory of the behaviour of materials under high pressure proposed by P.Savic 
and R. Kasanin.lt is based on the Coulomb interaction supplemented by a micro- 
scopic selection rule and a set of experimentally founded postulates. The theory 
has been applied to cases ranging from DAC experiments to calculations of 
models of planetary internal structure. 
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Introduction 



The study of materials under high pressure (and temperature) is impor- 
tant in a variety of situations in physics, astrophysics and related sciences. These 
range from highly exotic examples such as the early Universe, to laboratory ex- 
periments performed in diamond anvil cells (DACs). Determining theoretically 
phase diagrams of solids under high pressure is an extremely complex problem in 
statistical mechanics (for example [1] ).It involves the choice of a proper Hamil- 
tonian of the system under consideration, the calculation of the free energyand 
then the determination of the regions (or points) in the parameter space in 
which the thermo-dynamical potentials become non-analytical functions. 

The purpose of this paper is to present briefly the main physical ideas and 
examples of applicability of a particular semi-classical theory of the behaviour 
of materials under high pressure.lt was proposed by.P.Savic and R.Kasanin [2] 
and nicknamed the SK theory for short. 

Initiallythe idea for work on this problem dates back to a paper by Savic 
[3], which had the aim of exploring the origin of rotation of celestial bodies. It 
emerged from this work that rotation is closely related to the internal struc- 
ture, and that a theory of the behaviour of materials under high pressure was 
needed to explain it. The theory which he developed with Kasanin proposes 
a simple algorithm for the calculation of phase transition pressure in mate- 
rials subdued to pressure. Recent work [4] gives the possibility of establishing 
the equation of state and calculating the bulk modulus and its first pressure 
derivative. Apart from laboratory studies, the SK theory has found applications 
in astrophysics ( such as [5]- [9] and references given there). Note that astrophys- 
ical results close to those obtained within SK have been derived by Vasiliev [10] 
using different methods. 

This paper has two more sections: the next one contains an outline of the 
basic premises of the SK theory, while the third one is devoted to a brief presen- 
tation of applications of this theory to laboratory experiments. 



2 



Basic premises of the SK theory 



The object of study in the SK theory is a mole of any material under high 
external pressure. SK is based on the idea that increasing pressure leads to excita- 
tion and ionisation of atoms and molecules that make up the material. Translated 
into quantum mechanical terms,this means that increased pressure provokes the 
expansion of the radial part of the electronic wave functions in the atoms and 
molecules that make up the material. Such an idea might seem strange at first 
sight, because one may be inclined to think that high external pressure leads to 
a "crunch" .However, this problem has received a quantum mechanical treatment 
only about a decade ago [11] in a series of papers aiming at a theoretical expla- 
nation of the ruby scale. The mean inter-particle distance a is defined in the 
SK theory by the relation 



where Na is Avogadro's number, p the mass density and A the mean atomic 
mass of the material. One can now define the "accumulated" energy per electron 

as 



It can be shown that a as defined above is a multiple of the Wigner-Seitz 
radius [12]. 

Mathematicallythe basic premises of the SK theory are expressed in a series 
of 6 postulates which establish the ratio of densities and accumulated energies 
in successive phases [6]- [8]. Due to the fact that these postulates assume that 
jumps of the density occur in phase transitions, the applicability of the SK theory 
is limited to first order phase transitions. In all calculation within SK an integer 
index (denoted by i ) appears. Physicallyit has the meaning of the ordinal num- 
ber of a first order phase transition occuring in a material. The phase indexed by 
i = ends in the critical point. Without entering into details,the main practical 
result is the following expression for the phase transition pressure pt r 



N A {2afp = A 



(1) 



E = e 2 /a 



(2) 



Ptr = {, 



0.5101p*;i=l,3,5,... 
0.6785p*;i=2,4,6,.. 



(3) 
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where 



p* = 1.8077/3* (Vy 4/3 2 4t/3 Mbar (4) 

Pi = 3(a] /3 -l)/(l-l/a i ) (5) 

and 

_ r6/5;i=l,3,5,.. 

ai ~ l5/3;*=2,4,6,.. V°J 

The symbol V denotes the molar volume of the material under standard con- 
ditions, and i is an integer index numbering first order phase transitions which 
occur in a material. Physically realizable values of this index are determined by 
the selection rule 

E*+Ej = E* (7) 

where Ej is the ionisation or excitation potential and Eq and E* are pressure 

dependent characteristic energies of the specimen,which take into account only 
the Coulomb part of the inter-particle interaction potential. In another line of 
development, recent work [4] has given expressions for the equation of state and 
implicitly for the bulk modulus and its first pressure derivative within the SK 
theory.lt was shown there that the equation of state has the form: 
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where Na denotes Avogadro's number, A is the mean atomic mass of the 
material, p the mass density and other symbols denote various constants calcu- 
lable within SK. 
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Applications in laboratory experiments 



The algorithm proposed by SK for the calculation of the phase transition 
pressure was applied to 22 different materials [7], [8]. Apart from 20 materi- 
als for which values of phase transition pressure were known experimentally 
or from various theoretical calculations (CqH(,]CH4]CC14]Cqq]CsI]CF 2 = 
CF 2 ; CdS; Al; Ba; CaO; Kr 2 ;Ne; LiH; Te; Te0 2 ; Mg 2 Si0 4 ; Pb; S 2 ;AIP0 4 ; 

GdAlOs : Cr 3+ ) the SK algorithm was applied to hydrogen and helium. 
They were included in the calculation because of their astrophysical impor- 
tance, although the actual existence of phase transitions under high pressure 
in these two materials has not yet been confirmed in static high pressure ex- 
periments. However, metallic-like behaviour was obtained in shock-compression 
work at a pressure lower than the theoretically predicted value [13], [14]. The 
relative discrepancies between the values of the phase transition pressure cal- 
culated within SK and those existing in the literature vary between nearly 0% 
and 30%. These differences are due to a variety of factors. 

For example, the precision of experimental data and of the input parame- 
ters in the calculations account for approximately ±10%. An important source 
of the discrepancies is the form of the interparticle potential. The SK takes into 
account only the pure Coulomb part of the interparticle potential, without the 
contributions of the charge distribution overlap and of the dispersive and repul- 
sive forces. Now, the relative contribution of the dispersive and repulsive forces to 
the full inter-particle potential is minimal for C,H,N,O.Interestingly,the discrep- 
ancies between the SK and experimental values of the phase transition pressure 
is also minimal for the hydrocarbons [7]. 

Taking into account the existence of the charge distribution overlap gives 
rise to three additional terms in the expression for the accumulated energy.The 
sum of these terms is pressure dependent and it can be positive, negative or zero 
[7] .The existence of these terms induces an error in the calculated values of 
the phase transition pressure, and the magnitude of this error is also pressure 
dependent. Due to space limitations, we have here outlined the influence of just 
two factors which contribute to the the relative discrepancies between the SK 
values of the phase transition pressure in various materials and those obtained 
experimentally. A detailed analysis is available in [7]. 

Instead of a conclusion, several comments about the SK theory see to be in 
order. The big advantage of this theory is its physical and calculational simplic- 
ity.On the other hand, this simplicity necessarily induces discrepancies between 
the experimental values of the phase transition pressure and those calculated 
within SK. Pushing the reasoning further, the existence of these discrepancies 
opens up the possibilities for improving the basic assumptions of SK by render- 
ing them more complex and physically realistic. 
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One of the possible directions for improving the theory would be the in- 
clusion of non-Coulombic terms into the expression for the accumulated en- 
ergy.Some work along these lines has already started. 
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